arXiv:1504.07635v2 [hep-th] 23 Nov 2016 


Prepared for submission to JHEP 


Drag phenomena from holographic massive gravity 


Matteo Baggioli,'^ Daniel K. Brattan.^*^ 

Departament de Ftsica and IFAE, Universitat Autonoma de Barcelona, 

Bellaterra, 08193, Barcelona, Spain. 

^Physics Department, Teehnion- Israel Institute of Technology, Technion City - Haifa, 
32000, Israel. 

‘^Department of Mathematics-Physics-Computer Science, University of Haifa at Oranim, 
Qiryat Tivon, 36006, Israel. 


Abstract: We consider the motion of point particles in a strongly coupled field theory 
with broken translation invariance. We obtain the energy and momentum loss rates and 
drag coefficients for a class of such particles by solving for the motion of classical strings 
in holographic massive gravity. At low temperatures compared to the graviton mass the 
behaviour of the string is controlled by the appearance of an exotic ground state with non¬ 
zero entropy at zero temperature. Additionally we find an upper bound on the diffusion 
constant for a collection of these particles which is saturated when the mass of the graviton 
goes to zero. 


^E-mail address: mbaggioli@ifae.es 
^E-mail address: danny.brattan@gmail.com 



Contents 


1 The theory of drag 


3 


2 Holographic massive gravity 

2.1 The heat capacity 


4 

8 


3 Motion of a string in holographic massive gravity 

3.1 The low velocity dispersion relation 

3.2 Steady state motion at constant velocity 

3.3 Late time unforced motion 


10 

11 

14 

19 


4 Conclusions 


23 


A Analytic results for the mobility 


25 


B Analytic results for the decay constant 


26 


Holography is an extremely successful tool for computing physical observables in 
strongly coupled theories by considering weak perturbations of a gravity theory in one 
higher dimension [1-3]. Recently, increasing interest and efforts have been devoted towards 
modelling condensed matter systems in an attempt to learn about real world strongly 
coupled systems. 

To more closely approach real materials incorporating momentum dissipation mech¬ 
anisms in the gravity framework is essential. Without such mechanisms weakly coupled 
and holographic computations often lead to unrealistic results such as the divergence of the 
optical conductivity at zero frequency [4]. Several methods to introduce non-conservation 
of momentum have been proposed such as turning on a lattice distortion at the boundary 
(either through the modulation of the chemical potential or a scalar source) [5—7], the use 
of probe brane models [8, 9], helical lattices [10, 11], striped phases [12] and finally disorder 


[13]. 


An alternative method that has proven to be a very convenient way of incorporating 
momentum relaxation is the introduction of a graviton mass [14]. Following the original 
work many generalizations have been pursued with the common aim of having a complete 
description of strongly coupled metals, insulators and the transitions between them. For 
example, models with linearly sourced scalars were introduced in [15] and further developed 
in [16-18]. A slightly richer model than the linear scalar uses a dilaton field [19] while others 
have exploited global symmetries [20]. 

While the framework of massive gravity reproduces Drude peaks [21] and the metal- 
insulator transitions associated with broken translation invariance, the precise interpreta¬ 
tion of the deformation in terms of which operators we turn on in the boundary theory 
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is unclear. For example, in many condensed matter systems non-conservation of momen¬ 
tum can be associated with several phenomena such as the existence of phonons or the 
introduction of impurities. We shall find, by examining the heat capacity and the drag 
coefficient, that turning on a graviton mass cannot be uniquely identified with either of 
these condensed matter phenomena alone. Nonetheless we do see features generic to our 
large class of massive gravity models, in particular, a pronounced peak in the diffusion 
constant. 

To try and illuminate the role of graviton mass we return to one of the simplest models 
in condensed matter theory: the Drude model. This describes the motion of a probe particle 
interacting with a medium. Generally a moving probe particle loses energy and momentum 
to the surrounding material leading to an effective viscous drag. In a weakly coupled theory 
momentum is typically lost by two body collisions and the production of massless particles 
via bremsstrahlung. Examining the drag coefficient can generically give clues as to which 
mechanism is dominant, what types of collisions are occuring and the nature of the emitted 
massless particles. 

If the medium has broken translation invariance, either by the introduction of impu¬ 
rities or periodic potentials, we expect there to exist additional channels into which the 
probe particle can lose energy. It is important to note that the behaviour of phonons and 
impurities in a material are well understood in the condensed matter community. As such 
holographic studies investigating these phenomena are unlikely to yield new physics. 

In [22, 23] a probe string is used to calculate the drag on a charged particle due to 
thermal Yang-Mills. The string is a probe of the spacetime and experiences momentum 
loss to the background without significantly affecting the latter’s energy. When the applied 
force is due to an electric field this is the holographic analogue of the Drude model. One 
endpoint of the string enters the black hole horizon, which describes the thermal state of 
the boundary field theory and acts to absorb energy and momentum supplied to the other 
endpoint. The second endpoint of the string in [22, 23] is attached to a D7-brane. One can 
compute the diffusion constant for quark flavour [23], 


ttVXT 

where A is the’t Hooft coupling, by considering the motion of the free string endpoint after 
it is given some initial impulse. In this paper we shall compute the diffusion constants for 
point particles in the field theories dual to holographic massive gravity noting that these are 
distinct from the U{1) charge diffusion constants typically considered [10-14, 16, 17, 19, 20]. 

We shall perform string drag calculations in a particular class of massive gravity 
models. To study the effect of momentum loss one would ideally like to study the motion 
of a string in the fully back-reacted spacetime dual to a boundary theory with a lattice. 
Doing this however requires a complex numerical calculation for the metric which obscures 
some of the physics. We could instead consider a model like that discussed in [24] which 
uses a weak scalar source in an Einstein-Maxwell-scalar theory to produce a modulated 
lattice. As was noted however the effect of this weak lattice is to give the graviton a 
mass precisely in the manner of [14]. Since we are interested in the generic features of drag 
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with momentum loss we shall allow arbitrary potentials for the translation breaking scalars 
which encode generic massive gravity theories. Thus we are taking a bottom-up approach. 
No consistent UV completion of these frameworks is known so far and determining such a 
string embedding represents an important research programme. 

As the physics we are interested in concerns IR features of the theory, knowing the 
UV stringy completion of our model is not necessary. It is very likely that a generic string 
theory embedding would contain more degrees of freedom, such as a dilaton field. We would 
expect the dilaton, which could couple to our Stuckelberg sector, to have a radial running. 
Non-conformality of the boundary theory and its consequences for the drag mechanism, due 
to the presence of homogeneous and isotropic (in the boundary theory directions) dilaton, 
have already been investigated in [25]. Its effects in the context of holographic theories 
with momentum dissipation have also been determined [19, 20, 26, 27]. In particular, 
a strong dilaton can greatly modify the IR solutions of the system leading to new (often 
insulating and gapped) phases of matter. In order to separate the consequences of breaking 
translational invariance from those due to the running couplings, whose effect we already 
know, we have assumed that the background dilaton is vanishing. 

We begin the paper by discussing the theory of drag and its relevance to the Drude 
model of charge transport. Subsequently we give a brief description of one model of mas¬ 
sive gravity and discuss contributions to the heat capacity at constant volume. We then 
compute the drag coefficient and energy/momentum loss rates to the horizon for a probe 
string. At low temperatures with respect to the graviton mass we find that these observ¬ 
ables are controlled by the formation of an exotic ground state with non-zero entropy at 
zero temperature. Our results also indicate the existence of an upper bound on the diffu¬ 
sion constant of a collection of the probe particles. 


1 The theory of drag 

Consider a particle moving with spatial momentum p in a viscous medium under a driving 
force /. We model the motion of this particle using the drag equation 

| = -7P + / (1.1) 

where 7 is called the drag coefficient. We shall generally examine situations where the 
speed of the particle is low enough that its motion can be considered non-relativistic so 
that the dispersion relation takes the form 

where T is the temperature of the medium. In this case we define the spatial momentum 
to be p = MeffU and for steady state behaviour at constant velocity we hnd v = //(yMeg). 
We are particularly interested in determining the drag coefficient of this particle. 

Observing the stationary velocity the particle achieves under an applied force / allows 
us to determine the combination yMgff. Moreover if we assume the particle is charged and 
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the applied force is due to an electric field then ( 1 . 1 ) is exactly 
conductivity^ of the Drude model is well known, 

q'^ TIT 


O'DC 


Meflf 


the Drude model. The DC 


(1.4) 


where r = is called the relaxation time, q is the unit of fundamental charge and n 
is the density of charge carriers. To obtain Meg we then consider the late time unforced 
behaviour of the particle with some initial momentum p( 0 ) and solve the drag equation 
( 1 . 1 ) to find 


p{t) =p( 0 )exp(- 7 f) . 


(1.5) 


A measurement of the ratio dtv/v will then allow us to determine 7 and subsequently Mgff. 

Moreover, if the force in (1.1) contains a random element we can consider the Brownian 
motion of the particle through the medium. Under suitable assumptions on the nature of 
the statistical force it can be shown that the diffusion constant D is given by D = T/(yMefj) 
(which is called the Einstein relation). Once again if the particle has charge q and we place 
it in an external electric field E it will accelerate. At late times the velocity of that particle 
is (u)terminai = ■ Defining the mobility p by (u)terminai = we find 

d = , ( 1 . 6 ) 

which is called the Nernst-Einstein relation. We shall use this expression later on to 
determine the diffusion constant of probe particles moving in a strongly coupled field theory 
assuming that the statistical force supplied by fluctuations of the black hole horizon satisfies 
the conditions necessary for (1.6) to apply. See [28] for further discussion. 


2 Holographic massive gravity 


We use a generalised model of massive gravity to supply the background upon which our 
string moves. This is inspired by [17] where massive gravity is realized in a fully covariant 
formalism through the explicit introduction of Stueckelberg fields. This is convenient in 
order to have easy control on the number of degrees of freedom and possible inconsistency 
issues for which we refer to [29]. The minimal covariantization in the present case requires 
a set of scalar fields transforming under an internal Euclidean group of translations 
and rotations in field-space (for further details about this procedure see [30]). Our massive 
gravity action is then recovered by a truncation to 2 -derivative operators: 


5 'hmg = Mp j -y/^ 


R 


^ ^ - mV ( 


( 2 . 1 ) 


^In holographic theories with a t/(l) charge sector ctdc can be shown to consist of two terms [19]: 

r^DC ~ repair -\- (Jdissipation . (1.3) 


The hrst term is due to quantum pair creation and can be non-zero at zero charge density. It corresponds 
to the conformal part of the conductivity. The second contribution is due to dissipative mechanisms and 
its value in the Drude model studied in undergraduate condensed matter courses is given by (1.4). 
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Figure 1. Various plots related to the solution, (2.2) and (2.3), of the background equations of mo¬ 
tion. Left: The emblackening factor against radial position in the spacetime for the linear potential 
V{X) = V at various temperatures. From left to right we set T ^ 0.77,0.57,0.44,0.35,0.29 and 
ma = 1. The blue dashed lines represent the zero graviton mass solution at the given temperature 
while the solid red lines have non-zero graviton mass. Right: The entropy density for vanishing 
(V(0) = 0, dashed blue line) and linear potential {V{X) = X, red solid lines) against temperature 
for ma = 1, 2, 3,4, 5 (bottom solid line to top). Note that for non-zero graviton mass, represented 
here by the example of the linear potential, the entropy is non-vanishing at zero temperature. 


with V some polynomial. This theory admits solutions where the scalars take on linear 
expectation values = a S( x* with the Kronecker delta^. On these solutions the mass 
matrix for the graviton is identified as oc ■ Note that the Lagrangian for the 

scalars has some similarities with the effective Lagrangian describing phonons [31, 32]. 
It is necessary for V {X) to be monotonically increasing in X to avoid ghosts and we choose 
V(0) = 0 so that asymptotically AdS spaces are a solution when the scalars vanish. 

As the Einstein equations are easy to derive we shall not present them here. A class 
of metrics which solve these equations are of the form 


ds'^ = gMNdx^dx^ = ^ 


dz2 

W) 


- + dxl 


( 2 . 2 ) 


where f{z), for general potential V{X), reads: 

3‘ 


f{z) = 


1- I — 


r-ZH 


dw 


- {miYz'^ / V 




(2.3) 


with a being related to some coupling between the scalars. For monomial choices of V(X) 
one of the parameters, m or a, is redundant and can be reabsorbed via an appropiate 
redefinition. For polynomial V{X) both the parameters are independent and can play a 
different role. Henceforth we set (. = 1. 

Considering (2.3) it can be seen that at z = zh we have a coordinate singularity. 
Despite the known potential issues of singularities in the context of black holes geometries 
in massive gravity theories [33], our solution seems to be totally healthy and smooth. This 
is because in our scenario the Stueckelberg fields have non trivial profiles only in the 

^This choice leaves the Ward identity for energy conservation unbroken while breaking spatial translation 
invariance. 
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spatial directions as compared to the Lorentz invariant case discussed in [33]. Therefore 
the invariant object is regular at z = zh and we do not need to concern 

ourselves with the remedy to these problems described in [33]. Another quick, but non- 
comprehensive test, is to consider the standard curvature invariants R^, RmnR^^ and 
RmnpqR^^^^ which all remain finite at z = zh- 

As all the emblackening factors of (2.3) have a zero at z = zh, where there is no 
curvature singularity, our spacetime can be interpreted as thermal with a temperature 
given by: 


T = 


3 

dvrzH 


1 - [a^zl) 


(2.4) 


To obtain this temperature we follow the standard procedure of continuing to Euclidean 
time and imposing smoothness of the manifold at z = zh through suitable periodicity 
conditions. This ensures that correlation functions on the boundary share this property, 
the hallmark of a thermal field theory. 

We can see from (2.4) that it is possible to set T = 0 at non-zero zh and have extremal 
black holes. It is worth noting that the non-trivial scalars, in particular a^, enter the first 
law like a magnetic held [15] (or conversely a charge in four spacetime dimensions). As a 
consequence it is not surprising that one can reach an extremal limit by dialing the graviton 
mass appropriately as the scalar hair will support the black hole. Indeed assuming that 
the entropy density is given by the area of the black hole horizon we hnd 


s = 2 ttMp/z'^ , 


(2.5) 


which is non-zero at zero temperature if there is a non-trivial solution to (2.4). The 
entropy density is also displayed in hg. 1 for vanishing graviton mass and an example 
potential V{X) = X. Indeed, for all but a few special cases where the emblackening 
factor has higher order degeneracies at z = zh the extremal IR geometry turns out to 
be AdS 2 X IR^; typical of extremal black holes. It would be interesting to study such IR 
extremal geometries and their instabilities. 

With the above remarks about the thermodynamics in mind we note that the position 
of the black hole horizon can be given as a function of the temperature and graviton mass 
by solving 


T = 


^.■PZY^ 


1 ^ T/ /' 2 2 N 

1 - (a %) 


< 


dvrzH 


( 2 . 6 ) 


for Zh- This is bounded above by the zero graviton mass result^. For the square root 
{V{X) = A^/^) and linear {V{X) = X) potentials, upon solving (2.6) for zn, we find the 

®To see this assume that both T and zu must be positive (the possibility of having negative temperature 
T < 0 goes beyond the scope of this paper). Thus, when this equation has a solution it must be the 
case that the quantity in square brackets is positive. Moreover because V{X) is monotonically increasing 
between zero and (to avoid ghosts) and R(0) = 0 (so asymptoticlly AdS spacetimes are a solution) 

the same term is bounded above by one. It equals one only when the second term is zero (for example when 
the graviton mass is zero). Thus proving the result. 
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simple and exact results: 

3 


While there is no such simple result for higher power monomial potentials, generically we 
see that for larger powers of X in V{X) the leading correction to the zero graviton mass 
result is smaller. At large temperatures {ma^ <C we can show that the leading 

correction to the zero graviton mass horizon position is 


' (l + ^) \ y(A) = 

> y[X)-X 
I 




3 

4^ 





2n 


lim (X-’^V(X)) + ... 
x^o ^ ’ 


( 2 . 8 ) 


13 = , 


(2.9) 


where we have assumed the leading term in the potential at small X is proportional to 
X”. Thus if <C T then we can ignore the effect of graviton mass and our spacetime is 
thermal AdS. 

At low temperatures we face an obstruction to determining zn in terms of small T, 
namely the massive gravity spacetimes potentially exhibit zero temperature black holes 
with zh given by the smallest real solution to 

'' (“'4) = ^ ■ (2.10) 

These black holes are interesting and particularly worthy of independent study. We note 
that (2.10) always has a real solution if 14(0) = 0 and V(X) is monotonically increasing. 
Let be the smallest real solution to (2.10). For the square root and linear potentials 
this is given by 


(0) ^ 1 f 3 , F(X) = XV2 
/3 \ ^/3 , V(X) = X 


( 2 . 11 ) 


For any potential V{X) we find the horizon position as a function of small T is given by 


_ .,(0) 


Zh = 


H 


1 + 


1 " 1 -^ 7 ^^) 

3 


v"(4”’ 


T + 0‘^{T) 


( 2 . 12 ) 


The above asymptotic expansions for the horizon position at small and large temperature, 
(2.12) and (2.8) respectively, will be particularly useful for determining the behaviour of 
the drag coefficient and momentum loss rates at the relevant extremes. 

There are of course many possible potentials we could consider here. We will often 
consider V{X) = X^, the monomial potentials^, with N = 112, 1,2,3. The choices N = 

^We have found that the features of generic potentials can often be identified from the features of the 
monomial potentials. We shall see that the smallest power of X in V{X) is almost always dominant at 
large temperatures and small p. In the opposite regime it is often the largest power of X that dominates. 
The interest in non-linear combinations, such as V (X) = X -|- kX^ comes from the fact that they give rise 
to metal-insulator transitions [17]. 
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Figure 2. Various plots describing the behaviour of the heat capacity for different potentials. 
The potentials displayed are the square root {V{X) = \fX, solid red), linear {V{X) = V, dotted 
brown), quadratic {V{X) = X"^, dashed orange) and cubic {V{X) = X^, dot-dashed grey). Left: 
The change in the heat capacity from the zero graviton mass value against temperature for various 
potentials and P = 15. The solid line which is always increasing represents the square root potential 
while the solid line (grey) with the turning point represents the cubic potential. Right: The heat 
capacity at low temperatures (T = 10“^) as a function of P for various monomial potentials. 


1/2 and N = 1 are due to the existence of the exact and simple results given in (2.7) and 
(2.11). For larger > 4 we are forced to use numerics exclusively but will also use these 
techniques for N = 2, 3 as the analytic results are not concise. For intuitive purposes we 
will consider V{X) = X + kX^, with k some positive real number, where appropriate. 


2.1 The heat capacity 

For fixed potentials there are two thermodynamic state variables in our background - 
volume and temperature. We expect the thermodynamic potential to be extensive in 
volume and thus we can drop the former leaving only a single scale - the temperature - to 
be concerned with. The response of the system to thermodynamic forces is encoded in the 
susceptibilities and for us in particular the volumetric heat capacity at constant volume, 


_ 2{AT:fTMl 
dT 3zh 


1 - 'f (V(«24) - 2a24^'(«^4)) 


-1 


(2.13) 


where s is the entropy density of the system defined in (2.5). We note that for thermal 
AdS, i.e. V{X) = 0, the heat capacity at constant volume is 


_ 2MT!a/2 

<-^1^=0 ~ 32 • 


(2.14) 


In a typical material there can be several contributions to the heat capacity coming from 
distinct fields, for example, charged impurities® and phonon excitations. We now examine 
the behaviour of this quantity for various potentials to ensure the thermodynamic stability 
of our background. 

Consider the monomial potentials V{X) = X^. The heat capacities are displayed 
in fig. 2. We can see that for the zero graviton mass result, (2.14), as the temperature 

®Generically uncharged impurities do not contribute to the heat capacity. 
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goes to zero the heat capacity vanishes. This remains true for non-zero graviton mass 
(left-hand plot of fig. 2) and the heat capacity is positive and ~ T at low temperatures. 
Thus the system, governed by the exotic ground states with non-zero entropy density, is 
at least locally stable to thermodynamic fluctuations at low temperatures [34]. This does 
not guarantee of course that it is stable to other types of fluctuations nor that we have a 
global thermodynamic minimum. See [35] for a general discussion of the thermodynamics 
of black holes in massive gravity. 

From the left plot of fig. 2 we can see that for the linear potential there is something 
like a saturation temperature where the effect of the graviton mass is largest. It may be 
tempting to interpret this as a Debye temperature. However, we note that while calculating 
the difference between the heat capacity at zero and non-zero graviton mass certainly 
captures the change in this quantity due to spatial translation symmetry breaking it is 
difficult to argue that the additional heat capacity is due to the introduction of a new 
species, such as phonons, alone. 

Higher monomial and polynomial potentials have heat capacities that are reduced from 
the vanishing graviton mass value over some intermediate range of temperature (fig. 2 
bottom right and top left). Naively one would expect the addition of extra modes to 
increase heat capacity. This result suggests that some of the entities that contribute to 
the zero graviton mass heat capacity require more energy to excite as the graviton mass is 
increased. This would offset the introduction of more modes and lead to the reduced heat 
capacity observed. 

For square root and linear potentials the heat capacity can be expressed analytically 
as 


Cy 


Cu|m=0 ^ 


l+6( 


' AnT ’ 

air) 


V{X) = XV2 
V{X) = X 


(2.15) 


At large temperatures compared to /? we find the difference in the heat capacity from its 
zero graviton mass value (2.14) is 


( o \ 2n—3 

— j hm^ (A-’"H(X)) + ... (2.16) 


except when n = 1 where the next to subleading term must be considered. It is clear that 
for n < 1 at large temperature the change in the heat capacity increases with increasing 
temperature. When n > 1 however the opposite is true. This can clearly be seen in the 
left plot of hg. 2. This expression confirms for large temperatures that the correction to 
the zero graviton mass result is negative when n > 1 and positive when 1/2 < n < 1. 

Importantly, for the monomial potentials we consider and V{X) = X + kX^, with 
appropriate choices of k, the total heat capacity is positive and grows with temperature in 
the range of temperatures and j3 examined. Thus our background is at least locally stable 
to thermal fluctuations. 
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AdS boundary - z = 0 



Black hole horizon - z = zjj 


Figure 3. A representation of a string stretching from the horizon to the boundary and moving 
at velocity v. The boundary of AdS sits at z = 0 and the black hole horizon at z = zh- A string 
stretching to the hyperplane z = z„i from the horizon is given by truncating one that stretches all 
the way to the boundary. The string drags behind the “free” endpoint and falls into the black hole 
horizon at zjj. 


3 Motion of a string in holographic massive gravity 

Now that we have discussed the background we turn to calculating the motion of a string 
in it. We wish to consider a string dangling into the black hole horizon with the other end 
fixed to a constant z = Zm hyperplace (see fig. 3). We shall assume that the translation 
breaking scalars only affect the string through the pullback of the bulk metric to the 
worldsheet. Thus the string action is just 

5 = -r(o) I d^a ^-det igMN{X)daX^dbX^) , (3.1) 

where T^g) is the string tension, (a) are the embedding scalars and (t“ = (r, a) are the 
worldsheet coordinates. The equations of motion for the string take the form 


0 = drTT'Jy + daTT^ (3.2) 

where vr^^, vr^ are given by varying the Lagrangian density with respect to drX^ and 
dfjX^ respectively. Assuming the string action takes the form (3.1) may be inconsistent 
if we attempt to embed (2.1) in string theory but as an exercise this allows us to isolate 
how much the change in the background metric from thermal AdS affects the string and 
subsequently the drag on the point particle. 

We shall always choose static gauge a = z, t = t and let the only non-trivial string 
profile be x^{z,t). The Lagrangian density is 


£ 



1 + f{z){dzx^) 


(9iX^)2 

IW 


(3.3) 
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and the canonical momenta are 


/ 


TT' 


tt; 


V 


TTi 


vri 


vri 


T"2 

^(0) 

z^C 


^(0) 

z^C 


( 


\ 


vr!; 


dtx^{z,t)dzx^{z,t) 

/(z) 

1 + /(2;)(a^xi(2;,t))2 

dtx^{z,t) 

/(^) 

/ , _ {dtx^{z,t))^ 

fiz) 

-f{z)dtx^{z,t)dzx^{z,t) 
f{z)dzx^{z,t) 


\ 


\ 


The equations of motion for the string worldsheet embedding are 

f dtx^{z,t)' 

— Ot 


z^f{z)dz (J^dzX^{z,t) 


C 


= 0 


(3.4) 


(3.5) 


(3.6) 


with the outgoing boundary condition at z = zh and where the string is subject a force 
on the constant z = Zm surface. We note in [23] that this force is supplied by a constant 
electric field on a L)7-brane whose embedding terminates at z = z^ due to the finite quark 
mass. Here we shall assume the force can be supplied by some mechanism and remain 
agnostic about its origin. 


3.1 The low velocity dispersion relation 


Generically we shall choose the string endpoint not dangling into the horizon to move at 
small velocities. This allows us to extract analytically the form of the dispersion relation. 
In this section we shall argue that it takes the non-relativistic form independent of the 
string length. Firstly we calculate the static energy of the string. Choosing the string to 
sit at a constant value of x^{z,t) = Xstatic solves (3.19) and gives a total static energy 

£^static(T) = / dz TTl{z,t) =T^q) (- - . (3.7) 

Jzm \Zra Z}iJ 

The momentum of this static string is zero. While the expression for the static energy (3.7) 
in terms of the horizon position and string endpoint is the same as that found for thermal 
AdS [22, 23] we note that the dependence of F^static on temperature, through zh(/ 3,T), is 
different. 

Now consider a string imparted with some initial momentum and allowed to relax. 
At very late times the drag on the string will make it static. Prior to this point we can 
consider it as having some small velocity which deforms the string from the static profile. 
Such an embedding can be described by x^{z, t) = Xstatic+(^dx^{z, t), where e is a parameter 
of smallness, and the equation of motion for the fluctuation becomes 


z^f{z)dz 


(J-^dzSx^{z,t)^ 


d^6x^{z, t) = 0 . 


(3.8) 


Solving for the fluctuation and using (3.4), which we do in appendix A, it is possible to 
determine the dispersion relation at low velocities. For ease we repeat the result, (A.6), 
here 

E = F,tatic(r) + ^ 
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Figure 4. Plots related to the mobility for different potentials against temperature and /3. The 
potentials displayed are vanishing {V{X) = 0, dashed blue), square root {V{X) = ^/X, solid red), 
linear {V{X) = X, dotted brown), quadratic {V{X) = X'^, dashed orange), cubic {V{X) = X^, dot- 
dashed grey) and mixed (y{X) = X -t-O.dX^, solid black). Left: The mobility against temperature 
with /3 = 1. Note that all the potentials at non-zero graviton mass share a crossing point. Right: 
The mobility against /3 for T = \. The mobility for the mixed potential at small /3 is approximately 
equal to the linear potential while for larger /3 it resembles the cubic potential. 

with Ecstatic given by (3.7). The quantity defined above is the mobility of (1.6). To 
produce this result we have assumed that the velocity of the string endpoint is given by 
v{t) = 17(0) exp(—yt), at late times, with 7 the as yet undetermined drag coefficient. We 
shall compute 7 in section 3.3 and thus show that this assumption is valid at late times. 

It is important to note before moving on that the mobility of (3.9) is inversely propor¬ 
tional to the entropy density of the field theory. Given that the entropy density increases 
for non-zero potential (see right plot of fig. 1 for example) the decreased mobility of fig. 4, 
for non-vanishing potentials, is consistent with the interpretation of broken translation 
invariance in the boundary theory. The reason for this is that when we remove the re¬ 
quirement that microstates of the system must preserve spatial translation invariance the 
number of allowed states increases and there are more channels by which the particle can 
lose momentum. 

Although the graviton mass and potential V{X) enter in the definition of quantities 
like the entropy and energy there is a single thermodynamic parameter for the background 
- the temperature - once we fix the potential. Thus we should compare the mobilities in 
theories with the same temperature normalised by /3. 

Due to the bound on the temperature defined by (2.6), the mobility is bounded above 


by 



(3.10) 


The expression for zero graviton mass is equal to /Xbound(T), displayed in (3.10). For the 
square root and linear potentials, using (2.7), we find 



/• 


-2 


V{X) = 


(3.11) 


V{X) = X 
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These and more general results for the mobility with various choices of the free parameters 
are displayed against temperature in hg. 4. We see that the mobility generically decreases 
with increasing /3 and decreasing temperature. For all potentials we see that for low values 
of (3 with respect to T, on applying (2.8), 


Inbound (T) 3 




2n 


lim (X-^V(X)) + 

X-5>0 ^ ^ 


(3.12) 


where the leading power of X in V{X) at small X is assumed to be proportional to X'^. 
From (3.12) it is also clear that the higher the power of X in F (X) the weaker the correction 
to the case of zero graviton mass. We have observed that the situation switches at large (3 
with respect to T. Monotonicity of the potential with these two results is then sufficient 
to indicate the existence of a value of where any pair of mobilities will have the same 
value (see left plot of fig. 4). For the monomial potentials® this happens when zn = ^. It 
would be interesting to understand whether the existence of this point has some physical 
implications such as enhanced symmetries. 

From (1.6) and (3.9) we can additionally determine the diffusion constant for a col¬ 
lection of string endpoints which has the form 

D = r/r(r) = . (3.15) 

2 ( 0 ) 

Examples of the diffusion constant are displayed in 5. Except for the case of vanishing 
potential there is a peak in the diffusion constant at some value of temperature which 
moves to larger temperature as the power of the monomial increases. It is also interesting 
to notice that the diffusion constant goes to zero at T = 0 in the case of non-zero graviton 
mass where the exotic ground state is formed. 

Moreover, because the mobility is bounded above by the zero graviton mass result the 
diffusion constant is also bounded above by its zero graviton mass value. The existence of 
this upper bound is deeply related to the fact that V{X) is monotonically increasing (no 
ghosts) and that F(0) = 0 (asymptotically AdS spaces are a solution to the background 
equations of motion). A naive argument of what should happen to the diffusion constant in 
the presence of broken translation invariance would be that breaking a symmetry permits 

®The mobility can be used to partially classify the potentials V{X). Consider Vn{X) = W{X)^ where 
W{X) is a polynomial in X. The horizon position zh can be determined by solving 

2 

^IF Zii = AttTzh - 3 . (3.13) 

As we noted earlier, by the monotonicity property of the potential and assuming that at large P larger 
powers of X give greater corrections to the zero graviton mass result, there must be a crossing point for 
any two potentials. In particular consider Vn+i{X) and Vn{X). Because the right hand side of the above 
equation is independent of N it must be the case that 

W(q"4) = ^. (3.14) 

The crossing point is given by solving this equation and is independent of N. Thus all the potentials of the 
form Vn{X) = W{X)^ have a fixed value of the mobility at some value of zb_- 
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Figure 5. The diffusion constant for various monomial potentials against temperature with j3 = 0.7. 
The potentials displayed are vanishing graviton mass {V{X) = 0, dashed blue line with no turning 
point), the square root {V{X) = a/X, solid red), linear {V{X) = X, dotted brown), quadratic 
{V{X) = X^, dashed orange) and cubic {V{X) = X^, dot-dashed grey). Note that for the non-zero 
mass graviton as the temperature drops to zero the diffusion constant also goes to zero. This is 
distinct from the zero graviton mass behaviour. 

previously forbidden processes. We expect these to contibute positively to dissipation and 
thus the diffusion constant to increase. Clearly this cannot be correct alone and the results 
displayed in fig. 5 are consistent with the idea that there is an increased energy gap for 
processes that contribute to dissipation. This naturally makes the interpretation of our 
results in terms of new processes, like phonon production or interaction with impurities, 
more difficult because the change in the diffusion constant, drag coefficient and so on cannot 
be attributed to such a process alone. 

Were the reader to have picked up the paper at this point they may be a little 
surprised as the famous result for thermal transport in holographic systems places a lower 
bound on the diffusion constant. A bound for diffusion in strongly interacting systems 
with charge was recently proposed in [36] which our system violates badly at sufficiently 
small temperatures with respect to (3. See [37] for another model which possibly violates 
this bound. If we are to accept that the conjectured bound of [36] is indeed universal this 
perhaps suggests an instability in our system at low temperatures. To help decide which, 
if either, solution is correct it is necessary to better understand the interpretation of the 
string in the field theory dual to massive gravity spacetimes. 

We note in passing that many of the dynamical properties of a collection of particles 
are related to the mobility. It would be interesting to understand whether there are conse¬ 
quences of our bound for other transport quantities such as the thermo-electric transport 
coefficient recently investigated in [38, 39] for duals to massive gravity. 

3.2 Steady state motion at constant velocity 

For completeness we shall work at arbitrary velocity initially and only take the low velocity 
limit at the end of the calculation. We consider the following ansatz for the string profile 

x^{z,t) = vt + (3.16) 
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which describes a particle moving at constant velocity v in the direction on each constant 
z-plane. We substitute (3.16) into the equations of motion (3.19) and solve for ^{z) taking 
the string to stretch from the horizon to the boundary [z^ = 0). We must impose a force 
on the string equal to the momentum loss rate to achieve steady state behaviour. Shorter 
string profiles, where z^ > 0, can then be found by slicing the profile of the “full” string. 

The Lagrangian of (3.3), assuming (3.16), becomes 


L = -^Jl + f{z){d,iY- 


m ■ 


(3.17) 


From this expression it is clear that ^ is a cyclic coordinate so we define 


SL _ 1 f{z)dzi 


(3.18) 


where tt^ is a constant. This expression can be inverted to give a first order differential 
equation for 




fjz) - 

f{z) - 


(3.19) 


where the sign is fixed to — by requiring the string to trail behind the free endpoint. We 
note that because f{z) interpolates monotonically between zero and one and G [0,1] 
there is a value of ^ where f{z) — changes sign. If the string turns over at this point 
and returns to the boundary then this is not a problem. However the situation with two 
string endpoints terminating on the boundary describes a pair of interacting particles. 
For a single particle the string must past through the black hole horizon and thus the 
denominator must change sign at the same point (lest the induced metric of the string not 
be real). As such we must find solutions to the equation 


fiz^iv)) = 


(3.20) 


and set vr? = —t-tt- 

When the spacetime is thermal AdS it is straightforward to solve (3.20) to find 


Z*iv) = ZH (l - v"^) 


zh(1- 


(3.21) 


However there are corrections for the massive gravity spacetimes and solving (3.20) will 
generically require a numerical approach. In the case of small velocities we can however 
determine that 


and thus 


z*{v) = zn 


1 




7r| = ^ f 1 + + O^iv] 




ttTzh 


(3.22) 


(3.23) 
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It is sufficient to know dz^{z) from solving (3.20) to determine the energy and mo¬ 
mentum loss rates. Substituting our profile (3.16) into (3.18) we see that 9t7r* = 0. Thus 
TT^ is constant along the string following (3.2). Hence we can evaluate tt^ at any point on 
the string including at the horizon. Using (3.18) we can thus determine the loss rates to 
be 


vrf = -r(o)U7r^ , TTf = , 


(3.24) 


respectively. In the low velocity limit the momentum loss rate becomes 


_ 


'(0)^ 


1 + 


1 


27rr2:H 


+ 0\v) 


(3.25) 


It is important to now clarify that there is non-zero energy loss for the string even in 
the case of non-zero graviton mass due to the probe limit. The way that linear sources 
for our scalars were chosen in (2.1) means we are only breaking the subgroup of the full 
Lorentz symmetry responsible for momentum conservation. The Ward identity for the 
time component is kept untouched and energy is a completely conserved quantity in the 
background theory. 

For the square root and linear potentials we can determine the momentum loss rate 
analytically as /(z*(u)) = is at most quadratic in z^{v). Only the linear potential has a 
compact form however which we display for posterity 


z*{v) 


zh{T,/3) 


zh{T,(3) 

„ / 


47rT 


(1 - u^) - zu{T,/3f 
l-ZH(T,/3)3(f)' 
1 + 12 (^) 2-1 


21 1/3 


V 


6 


1^47rT J 


(3.26) 


(3.27) 


Noting that v enters as a square we see that only for quite relativistic velocities is there 
a signficant difference from the u = 0 result. We display the momentum loss for generic 
potentials in fig. 6. Beginning with the linear potential in the left plot we see that at large 
temperatures the momentum loss tends to the same large T behaviour i.e. T2. This 
large T behaviour is the same as the momentum loss of the zero graviton mass case, 


TT^ 


T>1 


vT^ 


(1 - u2)2/3 


(3.28) 


which is sensible because the effect of graviton mass is reduced by large temperature. At 
small temperatures however there is a non-vanishing momentum loss for non-zero graviton 
mass unlike in the case of zero graviton mass. Gravitationally this is because there are 
non-trivial solutions to (2.10) for the linear potential. In terms of the field theory we can 
see that the entropy density is non-zero at zero temperature. Thus there are states that 
the point particle can lose energy and momentum to. 
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Figure 6. Various plots of the momentum loss obtained from several potentials. The potentials 
displayed are vanishing (V(X) = 0, straight dashed blue lines), square root (V(X) = '/X, solid red), 
linear {V{X) = X, dotted brown), quadratic {V{X) = X^, dashed orange) and cubic {V{X) = X^, 
dot-dashed grey). Top left: The logarithm of the momentum loss rate for the linear potential 
against the logarithm of temperature at /3 = 1 and ln(r!) = —4.1,—3.1,—2.1,—1.1,—0.1 (bottom 
dotted line to top dotted line). The blue dashed line represents the zero graviton mass result. Top 
right: Logarithm of the momentum loss rate against the logarithm of temperature for various 
potentials with v = 0.1 and /3 = 2. The zero graviton mass line (blue dashed) has gradient two in 
agreement with (3.28) and we see at small temperature the loss rate tends to a constant for all other 
potentials. Bottom: The momentum loss rate against /3 for monomial potentials with ln(r) = —2 
and u = 0.1. It is clear that as the power of the monomial and f3 increase so does the loss rate. 


The other monomial potentials follow the same generic pattern. At large temperatures 
and small /3 we find that the momentum loss tends to (3.28). At the opposite extreme we 
find that the formation of the zero temperature black holes leads to a momentum loss 
approximately independent of temperature. Physically there is a non-zero density of states 
for the string to radiate momentum to at zero temperature as indicated by the non-zero 
entropy density. Again, we can understand this from the gravity side by solving for the 
position of the horizon at zero temperature 

AO) _ 

^ ■ /3 

As the power of the monomial potential increases the momentum loss at small temperatures 
increases (see top-left of fig. 6). Finally, from the physical point of view, increasing the 
power in the monomial potential should increase the strength of the interactions responsible 
for momentum dissipation and therefore increase the momentum loss itself. This can be 
seen from the bottom plot in fig. 6. 


(3.29) 
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Figure 7. The embedding of the string obtained from several potentials to be compared with 
fig. 3. The potentials displayed are vanishing {V{X) = 0, dashed blue line with no turning point), 
square root {V{X) = \fX^ solid red), linear {V{X) = X, dotted brown), quadratic {V{X) = X^, 
dashed orange) and cubic {V{X) = X^, dot-dashed grey). Top left: The embedding of the string 
for vanishing graviton mass with v = 0.1 and temperatures T = 0.1,0.2, 0.3,0.4, 0.5 (bottom to 
top respectively). As the string asymptotes to the horizon position at z = za the coordinate 
diverges. Top right: The string embedding for various monomial potentials with /3 = 1, T = 1/10 
and u = 1/10. For the embeddings to be distinguishable we have chosen /3 ^ T. 


On the condition that the velocity is small the momentum loss has the following 
generic behaviour 


7r_ 


vs 


2/N 2 
vm ' a 


(3.30) 


where again s is the entropy density. This is a lower bound for the momentum loss in 
the velocity as for non-zero v we see that z*(n) < zp- Unsurprisingly the momentum 
loss increases as the velocity increases. The approximate relation for the momentum loss 
(3.30) at small velocity also indicates that for m > 1 as increases the momentum loss 
increases while for m < 1 it decreases. This can be seen in fig. 6 where the zero temperature 
momentum loss increases with increasing N. 

Having computed the value of vrg the embedding equation (3.19) must generically be 
numerically integrated to determine ^{z). At small z we see from (3.19) that 

dz^iz) = —z^TTg-v/l — + ... (3.31) 


and thus we require that ^(0) is a constant. Subsequently we can integrate (3.19) to 
determine the string profile 


x^{z, t) = vt — 



(3.32) 


which we display in fig. 7. To describe a string stretching from a constant z-surface, z = z^, 
to the horizon that starts at = 0 at t = 0 we simply slice our embedding and use the 
coordinate shift 


/ ^iii 

/ ds 


f{s) - 


=0 f{s)y f{s) 


(3.33) 
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As we noted the canonical momenta along the string do not vanish. Following [ 22 , 23] a 
force must be applied to the string endpoint at z = Zm so that we cancel the momentum loss 
and subsequently satisfy the variational principle. Setting the force equal to the momentum 
loss rate and substituting into the force balance equation yields 


/ - ^p = 0 


f • V 

7-^eff = 

/J iZi 


'( 0 ) 


L 


1 + 


2ttTzi{ 


+ 0^{v) 


(3.34) 


where here • indicates the usual dot product and we have assumed low velocities so that 
p = Mesv. We see that the expression for to the lowest order in iP matches that given 
in (3.9). 


3.3 Late time unforced motion 


Earlier, in section 3.1, we described the behaviour of a particle at late times but neglected to 
determine the decay constant 7 as this required a numerical quasi-normal mode analysis. 
We now seek to remedy this omission. Using the Fourier decomposition of (A.l) and 
substituting into (3.8) leads to 


z^f{z)dz (J^^dzSx^{z,t)^ + uj‘^6x^{z,t) = 0 . 


(3.35) 


We will solve this equation numerically for the lowest quasi-normal mode as a function 
of the temperature and /3 by imposing outgoing conditions on the past horizon and then 
searching for io such that dz6x^{z,t) = 0 at some choice of Zm- There are also two regimes 
where we can find analytic results to compare against these numerics: small and large 
static energies. 

Before performing this analysis we note that our major results so far (namely the 
mobility and momentum loss rates) in sections 3.1 and 3.2 were independent of the string 
length. However the drag coefficient is dependent on this quantity. We have already argued 
we should compare theories of massive gravity at the same temperature. Now we make the 
choice to compare strings in those theories with the same static energy (3.7). For a choice 
of mass and potential picking the temperature amounts to fixing zu. Comparing strings 
of the same static energy then fixes Zm- Thus the free parameters are the temperature T, 
static energy Fstatic, /? and the potential. 

We display the results of numerical computations of the drag coefficient for various 
potentials in fig. 8 and fig. 9. We begin hrst with vanishing graviton mass summarising 
the results of [22, 23]. For small static energies the decay constant of the string is linear in 
the temperature and independent of Fgtatic, 

7L=o = 27rr 

For large static energies the decay constant is 

47rT rp2 

3 ^ 

1 J_( -^static I -^'static 

+ 4.T To) ) 



(3.36) 
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Figure 8. Various plots related to the drag coefficients for vanishing graviton mass {V{X) = 0 
dashed blue) and linear potential {V{X) = X, dotted brown). Top left: The logarithm of the drag 
coefficient at zero graviton mass against IniJstatic for InT = 0,1,2 (lowest to highest dashed blue 
lines). The numerical results are given by the dashed blue curves and tend to a constant at small 
i^static spaced linearly in InT. The black straight dashed lines are given by the analytic results 
in (3.37) and are thus separated vertically by 2InT. Top right: The drag coefficient against /3 
for the linear potential with InTgtatic = ~2, —1,..., 4, 5 (top to bottom line) and T = 1. We see 
that as Tstatic becomes small or large the drag coefficient becomes approximately independent of 
/3. Bottom: The difference between the drag coefficient for the linear potential and that at zero 
graviton mass against temperature. We fix Tstatic/r'(o) = 1 and choose /3 = 1,2, 3, 4, 5 (bottom to 
top line respectively). 


On fixing the temperature and varying the static energy we see in the top left plot of fig. 8 a 
smooth interpolation between a constant and 1/ii^static behaviour at small and large ^static 
respectively. 

Turning now to the linear potential, in the bottom plot of fig. 8 we display the drag 
coefficient for various values of the graviton mass, with ii^static/2~'(o) = 1, against temperature 
normalised by the zero graviton mass result. The difference is always positive indicating 
that there is an increased drag when the graviton mass is non-zero. As expected when the 
temperature increases this effect diminishes as the spacetime begins to resemble thermal 
AdS. When the temperature is small with respect to both /? and T^static we find that the 
drag coefficient tends to a constant. It then increases with increasing temperature before 
peaking and subsequently decaying in a manner consistent with the high temperature 
result. In the top right plot of fig. 8 we can see that as /3 increases the drag coefficient 
also increases. The amount that the drag coefficient grows by is dependent upon the static 
energy of the string with smaller energies and larger energies being less affected by the 
graviton mass. 
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Figure 9. Various plots related to the drag coefficients for the monomial and polynomial potentials. 
The potentials displayed are vanishing {V{X) = 0, dashed blue), the square root iV{X) = \/X, 
solid red), linear {V{X) = X, dotted brown), quadratic {V{X) = X^, dashed orange), cubic 
{V{X) = X^, dot-dashed grey) and mixed {V{X) = X -b solid black) Top left: The drag 
coefficient minus its zero graviton mass value as a function of F^static for various potentials with 
T = 1 and f3 = 1,2,3. Note that in both extreme limits of ifstatic the difference goes to zero 
in accordance with (3.38). The drag coefficient minus its zero graviton mass value as a function 
of temperature for various monomial potentials with ifstatic = 1 and /? = 1. Top right: The 
drag coefficient for the non-linear potential V{X) = X + kX^ with n/rri^ = 1, ifstatic = 1 and 
/3 = 1. We note that for low temperature the coefficient asymptotes the highest power behaviour 
of the potential; while for large temperature its behaviour is governed by the lowest power of the 
potential. This feature is quite universal and it is the inverse of the ^ dependence case where the 
highest power dominates in the large /3 limit. Bottom Left: The difference in the drag coefficient 
from it’s zero gravitom mass value against temperature for various potentials with ifstatic = 1 and 
/3 = 5. Bottom Right: A plot demonstrating the scaling behaviour of the drag coefficient in 
temperature for various monomial potentials with ifgtatic = 1 and /3 = 5. For large temperature 
the drag coefficient has a behaviour ~ jg reflected in the fact that the Vgcaiing ~ V at 

large temperatures. 


In the top right plot of fig. 9 we display the drag coefficient for various potentials 
against /3. Again we see that as /3 increases the drag coefficient also increases. When /3 
vanishes we get the zero graviton mass value for the drag coefficient which is generally 
dependent on the static energy and the temperature. For small values of /3 the larger the 
power of the monomial potential the larger the drag coefficient. This switches at some 
finite value of /3 so that potentials with greater powers have smaller drag coefficients. 

From the top left plot of fig. 9 we also note that at large or small Fistatic the dependence 
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of the drag coefficient on the potential diminishes just as was seen in the top right plot of 
fig. 8. In the bottom left plot of fig. 9 we another generic feature of the drag coefficient - a 
peak in the temperature. It is interesting that the existence of this peak is a generic finger¬ 
print of the model which does not rely strongly upon the choice of a particular potential. 
Additionally, for sufficiently large temperature, the difference in the drag coefficient scales 
as T~‘^^ (where N is the power of the monomial potential) as shown in bottom-right of 9. 

Having discussed our numerical results we would like to interpret them. For small 
static energies at non-zero temperature we use the analytic results of appendix B, namely 
(B.2), to show that graviton mass has no effect on the decay constant. This is consistent 
with the top right and top left plots in fig. 8 and with dy going to zero for small F^static 
in the top left plot of fig. 9. In terms of the gravity theory this is unsurprising because 
the near horizon geometry of (2.2) takes the same form for both thermal AdS and all our 
massive gravity spacetimes at non-zero temperature. The interpretation in the held theory 
is also relatively straightforward - namely the static energy is much less than the average 
thermal energy and thus the string is excited to the average thermal energy independent 
of the small differences in mass. 

For small and large static energies we use (B.7) to hnd 


7 = 


27rT , 
1 


.^'static ^ T 
.^'static T 


(3.38) 


Clearly in the limit of large temperatures compared to /3, where the graviton mass be¬ 
comes negligible, this reproduces (3.36) and (3.37). More generally, for Figtatic T, 
7 ~ s/Estutic ~ l/(li(F')£'static) where the behaviour of the mobility is displayed for various 
potentials in hg. 4. As we noted the mobility is bounded above at hxed temperature by 
the zero graviton mass result and is a decaying function as /? is increased and tempera¬ 
ture decreased. Similarly therefore the decay constant of strings with large static energies 
is bounded below at hxed temperature and generically increases with decreasing /3 and 
increasing temperature. 

For completeness we also determine the effective kinetic mass. In the limits of small 
and large static energy it has the analytic form: 


^ f 27rT4(/3,r) ’ ^static < T 

^(0) I F/gtaticj .^static T 


(3.39) 


This is compared to numerics in hg. 10. The left-hand plot shows the dependence of the 
effective kinetic mass upon the static energy for the linear potential V{X) = X at various 
temperatures and hxed (3 = 1. It is clear that as the static energy is increased we get the 
linear result given in (3.39). As the temperature is increased in the limit of small static 
energies the effect of (3 is diminished and we get the analytic result 


MeflF ^ 

^( 0 ) 9 


(3.40) 


This can be seen in the left hand plot of hg. 10 by the fact that for large temperature the 
brown lines at small Fistatic have approximately linear spacing. The right hand plot of hg. 10 
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Figure 10. Various plots for the effective kinetic mass. The potentials displayed are vanishing 
{V{X) = 0, straight dashed blue line), square root {V{X) = VX, solid red), linear [V{X) = X, 
dotted brown), quadratic {V{X) = dashed orange) and cubic {V{X) = X^ , dot-dashed grey) 
Left: The logarithm of the effective kinetic mass of the string against the logarithm of the static 
energy in the background of linear potential, V{X) = X, for (3 = 1 and InT = — 3, —2,..., 2, 3 
(bottom brown dotted line to top). At large and small Astatic with respect to temperature we see 
the result (3.39). Right: The effective mass a string against (3 for various potentials with T = 0.04. 

shows the effect of varying /3 for the monomial potentials at fixed T = 0.04 and Astatic = 1- 
Above some value of (3 the larger the power in the monomial, the larger the effective mass. 
Because we have chosen a relatively large Agtatic compared to the temperature as (3 is taken 
to zero the monomial potentials have a kinetic mass which is almost equal to the kinetic 
mass of the string at vanishing graviton mass. 

4 Conclusions 



Broken translation invariance 

Unbroken translation invariance 

V(X) - 0 

v(x) = Vx 

V{X) = X 

V(X) = n > 2 

c„/ Cylm = 0 

( 1 + iIt) 


T^ooi 

1 2(4,r)3T2 2 

V / 


Ai/Mbound(T') 

(1 + ^)-^ 

( + V 

1 <Ar 1 


^bound{^) = ( 4 ^) 

7 

(-^static ‘C T) 

- 

- 

27rT 


7 

(^static ^ 

- 

- 

rp2n 

~ - 

•^static 

4-7rT 

7U=o - X \ 

11 3 1 -^static 1 

+ '^'(0) ) 


Table 1. A representative sample of the results available in this paper. Cy is the volumetric heat 
capacity, fj, is the mobility, 7 the drag coefficient, T the temperature, T(q) the string tension. Astatic 
the static energy of the string, m the graviton mass, Mp the Planck mass and j3 a model dependent 
parameter which can be chosen freely. Results for more generic potentials can be found in the 
paper. 

We have calculated the momentum loss rates of probe strings in various solutions to 
a particular, very general, model of massive gravity. We have found that the mobility and 
subsequently the diffusion constant of the dual point particles is bounded above by the zero 
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graviton mass result. Moreover we have extracted the contribution of the graviton mass 
to the decay rate of a particle given some initial impulse. These results are summarised in 
table 1. 

We have argued that at low temperatures the behaviour of the string is governed by 
the existence of a set of non-trivial ground states. This is of course predicated on the idea 
that these ground states, which have non-zero entropy density at zero temperature, are 
stable. However, it is known that often such states are unstable. As these new ground 
states are quite exotic they are worthy of independent study which would simultaneously 
determine whether we can trust our results at low temperatures. 

We have seen that it is not possible to assign the additional momentum loss felt by 
a probe particle in field theories dual to massive gravity to a unique condensed matter 
phenomenon such as phonons or impurities. Our work thus indicates that extreme caution 
must be employed when using condensed matter terminology to describe the outcome of 
holographic experiments with massive gravity spacetimes. However, it has been argued [24] 
that a small graviton mass can be generated by a weak lattice in a gravitational version 
of the Higgs’ mechanism. Understanding the relationship between a probe string moving 
in a background with a non-weak lattice and the subsequent effective massive gravity 
theory should in principle help to clarify the role of different mechanisms, such as phonon 
production and interactions with impurities, responsible for the additional momentum loss 
observed in this paper. 

Finally, we made the simplifying assumption that our classical string was uncharged 
under the scalar fields. If the total action given by the sum of (2.1) and the string action 
is to preserve the shift symmetry —)• -t- so that we can treat the 4>^-helds as 

Goldstone bosons for broken translation invariance, the action of the string can have the 
more general form 

5 = -r(o) I d^a e/(^U9-*-)y_det (^gMN{X)daX ^, (4.1) 

where and f^\d^) are two unknown functions of ■ The function 

describes how the string is charged under the helds while r(o) exp(/^^)(cl<h)) describes 
some effective tension. It is also possible to add a Fradkin-Tseytlin term to the action^. 
Such a term would couple any non-trivial background dilaton to the worldsheet curvature 
and signihcantly complicate the equations of motion. As we have assumed that the back¬ 
ground dilaton vanishes we have not worried about such a term. Even in the presence of a 
non trivial dilaton prohle however two comments are in order: hrstly, the Fradkin-Tseytlin 
term would probably be an a' correction, which is certainly interesting but beyond the 
scope of the present paper. Secondly it is not clear from the holographic point of view 
whether such a term contributes to the on-shell action or whether it represents a pure 
counterterm®. That said, it would certainly be interesting to understand the effects of such 
a term. 

^We thank the anonymous reviewers for pointing out that such a term could have non-trivial effects. 

®We would like to thank E. Kiritsis for discussions on this point. 
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It may be interesting to ask how the physics changes when dependence on the scalar 
fields enters through more than the pullback of the bulk metric. We still expect to see non¬ 
zero momentum loss at zero temperature, if only due to the interaction of the string with 
the background scalars, although it is unclear to us whether our bound on the diffusion 
constant should survive. Moreover, it may be possible to motivate from a string model 
an effective held theory where the string is charged under the translation breaking scalars 
using methods similar to [24]. This would allow a precise interpretation of the string in 
terms of the boundary held theory as well as providing explicit forms for the undetermined 
functions in our string action. 
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A Analytic results for the mobility 

Consider the huctuation equation (3.8). Shifting to Fourier space with the conventions, 

(AT) 

we hnd the spatial momentum of the slow string is 
p{uj)= dz 7rl.{z-,uj) = e 

2 = 2ma,x 


/ d(jj 

— (5x^(2;; w) exp(zwt) 
27r 




lUJ zt, 


+ 0\e), 


(A.2) 


where we have regulated the momentum in the IR with an explicit cut-off Zmax which in 
principle we should take to be zr- Satisfying the outgoing condition at the past horizon 
requires /(zmax)d25x^(zmax; w) = za;(ix^(zmax; w) as Zmax -t ZH- Integrating p(u;) against 
frequency we see 


du! 


pit) = I p(w)exp(zwf) = ^(^(5x^(Zmax,t) 


27r 


(A.3) 


Turning now to the energy density of the string. 


4iz,t) = —^ 


2 / ( dtSx^(z,t)f 

2/(z) 2 


1 + " t))^ ) 0^(e) 


(A.4) 
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and integrating against the string length we find the total energy of the string is equal to 


E = ^static(r) + [5x^{Zma.y.,t)dt5x^{Zra^^,t)\ + 


2 r 2 

•^^max 


(A.5) 


where we have imposed the outgoing condition /(zmax)<9^5x^ = dtSx^ as —>■ -zh- 

If the free string endpoint is to satisfy the drag equation (1.1) it must be the case that 
p{t) = p(0) exp(—yt). Using the low velocity limit and integrating the velocity against time 
implies 6x^{t) = 5x^(0) exp(— 7 t) up to an additive constant we set to zero. Substituting 
into (A.5) then gives 


—Q, ^ 

E = E,t.tUT) + / + 0\p) , = yMeff = , (A. 6 ) 

2Mes{T) z^ 

with illstatic given by (3.7). We have yet to show that we can arrange for the endpoint to 
move according to p{t) = p( 0 ) exp(—yt), and determine the drag coefficient 7 , which we 
shall do this in section 3.3. 


B Analytic results for the decay constant 


We can extract some analytical behaviour for our drag coefficients in the regimes of small 
and large static energies. Beginning with the former we solve the fluctuation equation 
(3.35) in a near horizon expansion and find the derivative of the fluctuation is 


dzSx^ 


J_ / zn - 
zn V 


ibJ 

47rT 


\ r X .^ 9 / 



where 


Setting Zra 


32(7rT)^ + ZYif"{z}i) {2 ttT — iuj) 
?>2 {'kTYz}i {2t:T — ioj) 


zh imposes a constraint on the frequency and we find 


(B.l) 


uj = 2nTi . 


(B.2) 


This is the same result one finds for a short string in thermal AdS [ 22 , 23]. 

A second regime where we can find an analytic solution to compare with the numerics 
is the long string. To solve the fluctuation equation (3.35) assume that we can perform a 
small frequency expansion for non-trivial w, 

6x^{z-,uj) = 6xf'\z;co) + Ijj6x^^\z;oj) + oj'^6x^^\z;ij) + O^{co) . (B.3) 


We will find that non-trivial solutions exist only when a; is of order z^ with Zm ^ 1. We 
impose the outgoing condition on the past horizon and use an overall scaling of the linerised 
fluctuation equations 


0 = 


z^f{z)dz f^^(92(5x(o)(z;w) 


+ CO 




z^f{z)d:, (J^dJxl2)iz-,co)^ + 6 x\q^{z]uj) 


z^f{z)dz ( ^-^d^5x\j^{z;uj) 




(B.4) 
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to set 


lim 

ZS-ZH 



= 1. 


(B.5) 


Solving the equations (B.4) up to second order in w and calculating the derivative at 
we find (logarithmic term is not quite right) 


djx^(z;u;)l = u 


4 




lim 


L=Zm W^f{w) 


+ 


AttTz^ 


In 


ZH- Z 
ZH 


+0^iuj) . 


(B.6) 


It is clear that for a non-trivial solution we cannot truncate to order a; but must cancel 
the order a; and terms. Expanding the Neumann condition to leading order in small 
zh — Zra we See find 


iLi = iz^lz\. (B.7) 

This is small, as required by our expansion (B.3), on the condition that the string is long 
i.e. Zm 1. 

We see that the long string has fixed dependence on z^ and zh as long as the line 
element has the form (2.2). The dependence on Zm and zu is the same as that found by 
[22, 23]. However, once again, because zr depends on the temperature and graviton mass 
the drag coefficient of long strings is distinct between different theories of massive gravity 
at the same temperature. 
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